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Z.NIAC  COMPUTATION  OF  TWO-DIMENSIONAL  SUPERSONIC  NOZZLES 


Nathan  Gerber 


SYISBOLS 


X,  y 

q 

U,  V 

q,  9 

% 

a 

y 


M 


■K 

P 


2 2 
H = a -u 


K = -uv 

X,  |i 


P 


s 

A e 


rectangular  position  coordinates 
velocity  vector 

rectangular  velocity  components 

polar  coordinates  of  velocity 

speed  of  radial  flow  on  initial  circular  arc 

local  speed  of  sound 

ratio  of  specific  heats  (=1.U  for  air) 

Mach  number  (*q/a) 

inclination  of  bounding  ray  of  radial  flow 
curvature  of  nozzle 

characteristic  coordinates  (a  varies  on  the  patching 
characteristics ) 

L = a2-v2 
R = a(q2-a2)^/2 

direction  cosines  of  straight  line  (oc  = const.)  simple 
wave  characteristic 

length  of  0<  - const,  simple  wave  characteristic 
density 

x,  y along  characteristic  which  patches  the  transition 
and  simple  wave  regions 

arc  length  along  transition  section  of  nozzle 
parameter  governing  rate  of  growth  of  curvature  of 
transition  section 
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ABSTRACT 

A procedure  based  on  the  method  of  characteristics  is  described  in 
detail  for  computing  two-dimensional  supersonic  nozzles  with  continuous 
curvature  everywhere.  The  results  of  ENIAC  calculations  of  a series  of 
nozzles  by  this  method  are  summarized  briefly  . 
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SECTION  I 


INTRODUCTION 

This  report  will  describe  a nethod  of  determining  two-dimensional 
supersonic  nozzles  on  the  basis  of  piano  isentropio  flow  of  a perfect 
gas.  The  computational  procedure  is  based  on  the  method  of  character- 
istics (ref.  [i]  )• 


j 


The  procedure  commonly  used  for  designing  these  nozzles  (ref.  M) 
is  to  patch  a !' simple  wave"  onto  an  assumed  initial  radial  flow  to  obtain 
the  "cancellation"  region  which  leads  to  uniform  flow.  These  several 
regions  are  indicated  in  Fig  1.1,  which  shows  a representative  two-dimen- 
sional nozzle.  (Because  of  symmetry  we  need  only  refer  to  the  upper  half 
ox  the  diagram.)  The  radial  flow  originates  from  a two-dimensional  source 
at  point  0,  the  ray  OB  being  tangent  to  the  contour  TBJ  at  the  "inflection" 
point  B.  BD  and  DJ  are  characteristics  (or  Mach  lines)*.  The  problem  is 
the  determination  of  BJ • 

Nozzles  designed  in  the  above  manner  will,  in  general,  have  discon- 
tinous  wall  curvature  at  points  B and  J (Fjg  1.1).  While  this  condition 
is  permissible  for  fixed  throat  tunnels,  it  is  unsatisfactory  for  flexible 
throat  tunnels,  the  walls  of  which  must  have  continuous  curvature  every- 
where. The  process  used  in  the  calculations  described  in  this  report  avoids 


* See,  e.g.,  pp.  259  - 286  of  ref.  [3J  for  discussion  of  simple  waves  and 
Mach  lines. 


this  trouble  by  inserting  between  the  radial  f]  ow  and  the  simple 


J 


Fig  1,2 

wave  a transition  section  (EDFE  in  Fig  1.2)  the  wall  curvature  of  which 
varies  continuously  from  zero  to  the  value  at  the  beginning  of  the  can- 
cellation section.  The  transition  section  can  be  fitted  in  many  ways. 

The  particular  method  used  here  is  described  in  Section  II  by  eqn.  (2.6) 
and  its  accompanying  discussion, 

A series  of  calculations  described  in  Section  IV  were  performed  on 
the  ENIAC  to  provide  design  data  for  the  13"  x 15>"  flexible  nozzle  wind 
tunnel  (No,  1)  now  under  construction  for  the  Supersonic  Wind  Tunnels 
Branch  of  the  Exterior  Ballistics  Laboratory  * These  shapes  were  selected 
from  a three -parameter  ( qQ,  9 , ^0)  family  of  nozzles.  The  particular 
choices  of  parameters  were  guided  by  factors  such  as  physical  limitations 
on  tunnel  lengths  and  permissible  plate  curvatures.  The  object  was  to 
get  a continuously  varying  family  of  nozzles  from  which  could  be  obtained 
reasonable  shapes  of  curves  ("jack  displacement"  vs.  Mach  number)  for  the 
design  of  cams  governing  automatic  settings  of  the  jacks  which  hold  the 
wall  in  place . 

To  obtain  the  final  family  of  eighteen  nozzles  required  about  one 
hundred  and  ten  trials  at  about  forty-five  minutes  per  case.  Clearly 
computations  on  this  scale  would  never  have  been  undertaken  by  hand. 
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SECTION  II 


EQUATIONS  AND  BOUNDARY  CONDITIONS 

The  equations  describing  the  two-dimensional  flows  are  (ref.  [l],  P.12) 

K Bu/Bx  + K (ou/By  + Bv/Bx)  + L Bv/cy  = 0 ^ 

Bv/3x  - 3u/3y  = 0 J 

2 2 2 2 
where  H x a -u  , K a -uv,  L s a -v  . 

The  flow  field  is  divided  into  four  regions  (Fig  2.1):  I — radial 

flew,  II  — general  two-dimensional  flow,  III  — simple  wave,  and  IV  — 
uniform  flow.  These  regions  are  separated  by  the  characteristic  curves 


y 


Fig.  2.1 

BD,  EF,  and  FJ:  BB'  is  a circular  arc. 

The  boundary  conditions  are 

a)  v*0  when  y*0j 

that  is,  the  air  flows  along  the  center  line  of  a cross-section  of  the 
nozzle . 


b)  v/u  = ff  (x)  along  a contour  f (x),  which  is  to  be 
determined  as  part  of  the  solution  to  the  problem.  This  contour  is  sub- 
ject to  the  following  conditions:  f'  (x)  = tan  0O,  f”  (x)  = 0 at  point 

Bj  f*  (x)  * 0 at  point  J;  f (x)  has  continuous  curvature-^  (x)  every- 
where between  B and  J . 
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We  first  establish  the  coordinate  system.  The  origin  is  deter- 
mined so  that  Xp  = 0,  y^  = 1 (Fig  2.2). 


D 


Fig . 2 .2 


The  initial  information  at  our  disposal  is  the  radial  flow,  with 
speed  qQ  on  HB' . Since  the  streamlines  coincide  with  the  rays  from 
point  0,  0 measures  the  inclination  of  both.  The  angle  Q Q is  divided 
into  n^  intervals  of  magnitute  cT  0 = 0 Q/n^,  where  is  an  integer. 


(2.2) 


Clearly,  on  equally  spaced  points  on  BB* 

x(n)  » esc  0Q  cos  (0Q  - n0Q/n1)  - cot  9Q 

y(n)  - esc  0Q  sin  (0Q  - r& 

u(n)  = qQ  cos  (9q  - 

v(n)  = qQ  sin  (90  - r^/i^) 

where  n * 0,  1,  2, ...,n  . Accordingly,  n = 0 for  point  B and  increases 
along  the  curve  up  to  rrj-  at  point  B'  . 


Vre  designate  c*  as  the  variable  along  HD,  the  Mach  line  patching 
together  the  regions  I and  II.  In  the  solution  of  problems  by  charac- 
teristics the  choice  of  characteristic  parameter  along  certain  bounding 
Mach  lines  is  arbitrary  (see  ref.  [1]  , p.  21).  Thus  we  are  able  to 
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choosea  ®n  along  BD;  that  is,  0(  will  vary  conti nuousDy  with  6 from  zero 
at  B to  n^=c<^  at  D,  and  at  any  point  P (Fig.  2.2),  where  9 =6  c(l-n/n^),CX 

will  be  equal  to  n. 

Since  the  holograph  of  the  plane  characteristic  BD  is  a Prandtl- 
Busemann  epicycloid,  the  velocity  distribution  on  BD  is  described  by 


9 - 9 » f(M)  - f(M ),  where 


f(M)  a cos"1  (1/M)  - ]/(r+l)/(r-D  tan"1  (/l^-l  /(T^l)/(r+l)), 

r 1 / o / r o?j  (2 «3) 

(see  ref.  [_!;],  p.  26)  and  M(q)  =\j2q  / [0''-l)(l-q‘‘)J  > also  Mq  =>  M(qQ)  . 

To  determine  the  curve  BD,  consider  the  streamline  OP  (Fig.  2.2) 

/ 

oF/$P  = yp/VP» 


UP 


= esc  9 


Thus  t ^P  = (yp/  ypt ) esc  9 


P and  P'  may  be  taken  as  two  points  on  a two-dimensional  nozzle  sur- 
face. Then  (ref.  [5],  p.  3U) 

yp/yp<  = (Ap/A*)  / (Ap'/A*) 


A/A-k-  is  the  ratio  of  the  cross-sectional  area  at  ary  point  in  a two-dim- 
ensional nozzle  to  that  at  the  throat  of  the  nozzle  (according  to  one- 
dimensional theory,  M-i  at  the  throat)  and  is  a function  of  q. 


A/A*  = (1/M) 


£+(*  -DM2}/  (>+!)] 


{y  +i)/(2  y -2) 


Since  qp,  («qQ)  and  qp  (cx)  arc  known,  fp  (or)  can  be  determined. 


Let  9 = 0Q(1  - n/n^)  = 0 0(l  /o<  «=  9 (or)  at  point  P. 


Then  ^ 

x(« of)  = fp  cos  (1  -cx/cXp)  90  - cot  9q 

y(CX)  = l p sin  (1  -«AX)  ®0 

u(crt)  = q(Of)  cos  (l  -CX/x^)  90 

v(<x)  » q(cx)  sin  (1  - ) 9Q 


(2.U) 
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Consequently  x(or),  y(C*),  u(cx),  v(Of)  determined  at  cx^+l  points  on  BD, 


We  transform  to  the  characteristic  (cj,  p)  plans  to  determine  the  flow 
in  region  II.  The  characteristic  equations  are  (i‘of.  [l]  )s 

"h 

(j'cx  =vy/*&,  yp  =cy/6p,  etc.) 


(K-R)  y<*  - Lxa=  0 
H yp  - (K-R)xp  « 0 
H Uo<+  ( K-R ) Vex  = 0 


(K-R)ud  + Lv 


P 


J 


(2.5) 


where  F - HL . Fig . 2 .3  shows  the  map  of  region  II  iti  the  character- 

istic plane . The  two  bounding  streamlines  BE  and  DF  map  into  the  lines  p“Of 


P 


Fig.  2.3 

and  respectively  (dp/d«=  1 on  them). 

We  arbitrarily  demand  that  the  inclination  of  the  nozzle  on  BE  be 
of  the  form 

9 * 9o  ^cx2  Ae  (2.6) 

where  A0>O  is  a prescribed  constant.  The  curvature  on  BE  is 
= d9/ds  = - 2<*  Ae/s^  • 

At  B,  of  a 0,  so  -f(  - 0.  Thus  the  condition  of  zero  curvature  of  the  nozzle 

at  point  B is  satisfied.  More  generally  9 = 0 -Of  g(cx),  where  g(Or)  is  an 
arbitrary  function  of  CK  , satisfies  this  condition. 

Since 

ds/dO(=  (Sx/SOf)  \/l  + (dy/dx)2  - (9x/9<x)  /l  + tan2  (0  -CX2A9), 
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the  curvature  is 

^ * - 2£7f  AeyY(dx/dO0  |/l  + tan2  (fl0-<*2  *^9)J  (2.7) 

The  flow  in  region  III.  is  a simple  wave;  that  is,  the  velocity  is 
constant  along  the  characteristics  Of  = constant,  which  are  straight  lines , 
Consequently  the  points  on  EJ,  the  remainder  of  the  nozale,  are  simply 
functions  of  Of  , In  particular,  0.  ■' 0 ) ■*  o,  and  is  the  speed  of 

the  final  uniform  flow,  1 


Fig.  2. I* 

We  wish  to  locate  a point,  say  T,  on  EJ  (Fig,  2,1*)  which  corresponds 
to  the  point  G on  EF.  Let  x =£  , y -Y)  on  EF.  Then  x^,  a 6 (°0 , yG  = Y)  (<*) 

The  slope  of  GT  is 

m(«)  = [(9y/ap)  / Ox/5p)  jG  ^ (K-P.)  / H (from  eqn.  (2,5)  ). 

The  unit  vector  tangent  tc  GT  is  ( X,n),  where 

X « - H/i/?  + (K-H)2  , \x  = - (K-R)  /h2  + (K-R)2. 

Then  on  EJ 

x(a)  = i (Of)  + X (Of)  r (Of)  - Xp 
y(0()  - IT)  ( PO  + n («X)  r (Of)  = yT 
u(a)  = = u^ 


(2,8) 


r (Of)  = GT 


v(Of) 


= Vr, 


(2.9) 


We  find  r(°f)  essentially  from  the  equation  of  continuity  by  equating 
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the  mass  per  second  of  air  flowing  across  GT  to  the  mass  per  second  across 
EG. 

“ mass/sec.  across  GT  = p(or)  |ja(o()  M-  (or)  - v(oc)  X(cx)J  r(cx) . 

(uji  - vX)  is  the  scalar  product  of  the  vector  q = (u,v)  and  the  vector 
(p,,  - X),  a unit  vector  normal  to  the  line  GT.  It  is  also  the  component 
of  the  velocity  normal,  to  the  characteristic  GT,  which  is  equal  to  a (or)  f 
the  speed  of  sound. 

From  ref.  [l]  , p.  12,  p = p (1  - q2 (p  being  constant 

o 2 ^9  % 

for  the  flows  considered  here;  h*l),  and  a = (V'-l)(l-q")/2 . Hence 


QgS  - P0  (1  - q2 )1/(^~1)  a (ex)  r (or).  

A unit  vector,  normal  to  the  curve  EF  is  {~  H / f \ V, 

j a r\vc  r » 

oc/l/^f  + );  and  QgG  s mass/sec.  across  EG  « V | pO^-  u 


*W  )£&  + y)w£ 

2l 

V ^»OC  -r  l/OT  j uo  . 

But  ds/dOf  =*  4*  h)cx2  ; 


Oof)/ 


so 


QEG  = po  uy)or)  d<*. 

E 

Equating  QgG  and  C^.-p,  we  get  (taking)'®  l.i|) 

r (a)  - \j£  (1-q2)2’5  (v^-  «^)  d«j  /[a(l-q2)2  *5] 

According  to  eqn.  (2.0)  ^LAK-RjJ^  = (K+R)£CT  /H  ; so 

pot 

^«E  (l-q2)2,5[v  - (K+R)u/‘H]E«dor|  / ja(°0  fl-q2(0()]2*y 

j-  J ' -}  (2.10) 

Along  EJ  dy/dx  = v(o<)  /u(oi),  dV/dx2  = |d(v/u)/d°'j  /(dx/do). 


finally 
r (or) 


The  curvature  is 


-fc  jVy/cbr2  j J |l+ (dy/dx) 2 J ^ = (u/q)3  d(v/u)/d«]  /(dx/dor) . 

At  point  E,  where  r = 0,  we  note  that  (by  eqn.  (2.9)) 
(djc/dcx)]E  = (d^/doO]  E + X (dr/da)]  £ 


(2.11) 

(2.12) 
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p 


By  eqn.  (2  .1C)) 
(1 


a dr/d<J  e’(I-‘j2)1/P  /•*=•_]  E -udr)  /d0]  E} 


This  leads  to  the  relation 


dr/dofj  E 


q2X  / (a2  -u2)  ^ d £ /dCX 


which  upon  substitution  into  eqn.  (2.12)  loads  to 

- ' ^2\u  \/q  -a 2 / (a2  - u2)Jd^  / da 


dx/d<X 


E 


By  eqn.  (2.5),  X du/'da  + p dv/da  = 0,  and 
2 / "o  ^ 2 

d(v/u)  / do(  = (uvor  - vu  )/u  = - ( v q'  - a / jjiu  ) uc< 


(2  ..13) 


Substituting  this  relation  and  eqn.  (2.13)  into  eqn.  (2.11)  we  finally 
obtain 


u2)  /(2q^ii) 


(u0f  / )g 


(2.1U) 


Actually,  the  length  of  the  transition  section  is  not  known  a priori. 
Tn  the  computation  at  each  point  on  BJ  (Fig.  2.1)  both  ^ (eqn.  (2.?))  and 
(eqn.  (2.14))  are  computed.  When  a point  E where  ^ is  foun(i, 

generally  bv  interpolation,  the  characteristic  EF  is  constructed  and  the 
cancellation  section  thus  started.  To  achieve  a desired  design  Mach  number 
at  the  uniform  flew  region  several  trials  must  generally  be  made  with 
various  initial  qQ*s  on  BB1 . 
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SECTION  III 


COMPUTATIONAL  PROCEDURES 

The  initial  data  for  determining  the  flow  on  BD  (Fig.  2.1)  are 
available  on  IBM  output  cards  as  a result  of  a previous  ENIAC  compu- 
tation of  Frandtl -Meyer  flow.  (See  eqn.  (2  .3)0  A /A*  and  qQ  have  been 
tabulated  at  .0$  degree  intervals  of  = 9 + const.  The  input  deck  for 
a computation  is  extracted  from  this  main  deck  by  selecting  first  the 
card  with  q *»  qQ  (on  which  to  is,  say,  ajQ),  then  those  cards  with 
co  = co+cTa.co  = co  + 2<T  9 , . . to  = to,  + n,cf  9 = to  + 9 . Eqns  . (2  .U) 

then  yield  x,y,u,v  at  Oj+1  equally  spaced  joints  on  BD  (|3=0)  in  the 
characteristic  plane  ( Fig.  2.3). 


For  region  II  there  are  three  computational  routines:  the  stream 

process  for  points  on 


e 


Fig.  3.1 

DF  (Fig.  3.l)j  the  nozzle  process  for  points  on  BE,  and  the  general 
process  for  all  other  points*  The  second  order  iterative  method 
described  in  detail  in  ref.  [l]  is  used  to  solve  the  differential 


1U 


equations  numerically  . The  partial  derivatives  are  replaced  by  differences, 
so  that  for  a point  F (Fig.  3.2),  with  information  known 

0 

4 L M V 

« 4 • 

Q.  N 

• • 

U 


Fig.  3.2 

at  points  L and  U,  eqns . (2.£)  become 

(K-R)^(yp“yL)  " l$m('xp_xiP 

" (K“R)^(Xp-Xu) 

H$^ap“'V  + ^K“R^^Vp"Vt^ 

(K“R)$j(up“utP  + vw 

The  notation  means  that  H is  evaluated  first  at  Q,  then  in  the 
second  iteration  at  H.  HQ=(Hjj+Hp)/2j  HM«(Hp+Hp)/2  j Hjj=(Hp+Hy) /2,  where 
Hp  is  computed  from  up,  vp  determined  in  the  first  iteration.  This 
notation  applies  similarly  to  all  other  coefficients* 


0 
= 0 
0 
0 


(3.1) 


The  solution  (general  process)  to  eqns.  (3.1)  i: 
xp=(K-R)^m|h^]  yy  - Ov-R)^  Xu"|/d  - 

H#  {(1Wi)i*  yL  - l#*l}/d 
yp  " l#m  {h$n  7u  “ <K-RV 

(K“R)$N  {(K-R)j^  yL  - Xl}/D 

up  = l#n  \ + 'K"r^'Ivl}/^  " 
^K-R'£m  \}K~R)0  uu  + l0n  vu}/^ 
vp  = h(^m  {(k“rV  uu  + L^vu}  I ^ ~ 

(K“R)j^j  uL  + ('K"R^<^1Vl}'  / ^ 


V 


(3.2) 


It  should  be  remarked  that  the  ENJ.AC  coding  used  was  actually  adapted 
(to  save  coding  time)  from  existing  axi symmetric  flow  coding  by  suppressing 
a2v/y  terms . 
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where 


® = ^01  ~ an(^ 
h$m  L$i  " (k_r)$j  VIV"n;^l 

Similarly,  the  stream  process  computation  (DF  in  Fig.  3.1)  for 
point  F in  Fig,  3«3  is 


fi 


■>-Of 


Fig.  3.3 

y"f-  ~ yl_  ^ |"(K-K/ /l}  yL  ^ 

'*p*“l‘  {<>'-R)/h]  p'-i, 

The  equations  used  for  the  ncsalo  process  (BE  in  Fig,  3,1)  are 

Hy„  - (K-R)  x0  » 0 
P r 

(K-R)up  + L Vp  - 0 ^ 

' ’ dy/dx  * tan  (0— c x^A0) 

v **  u tan  (B.  -oi^AO) 

o j 


>-  cx 


Fig.  3Jl 


(3.3) 


(3.U) 
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i.a>J 


(3.5) 


xn  difference  form  they  are  (see  Fig.  3.U) 

[h/(K-R)1  (yp-ytJ)  - (xp-x,,)  = 0 

(■'ip  - ’^j)  + |l/(K-P.)J-  ^y(vp"v-j)  e 0 

(yp  - yA)  /(xp-xA)  - tan  M (9q  -of2A0)  - 0 

up  tan  (0O  -ot|a9)  - vp  = 0 _ 

where  taiVj  (»o  - cX 2 AO  ) = (1/2)  [jtan(0  of2A9 ) + tan  (9o  -Cx2A9)J  ■ 

The  solution  to  eqns.  (3.5),  which  locates  the  nozzle  points,  is 
Xp  = (1/A){[.H/(K-R)j  tj'q  yy  “ Xyj-  + 

(i/A)f[H/(K-R)]  j,q  [-rA.xA  tan„<eo  -or2Ae)]j 

yp  = (1/A)  ^-yA  + xA  tanM(9Q  -of2 A0 ) - - 

(1/A)  f-  [H/(K-R)]  y^xjl  tan^(9  O-O(2A0 ) 

”r  - <i/fc){[L/(K-a)]  V'Hi) 

Vp  • (l./B)  f[L/(K-K,,,,.  v„+uIr\  tan  (0  -M2A9) 

i i ww  u v I u r 


>■ 


(3.6) 


A 

B 


r ^ p 

|H/(K~R)|^q  tan  M (90  -cfA*  ) - 1 


- |l/(k- 


■)  P 

■R)j  tan  (9  0 -cXp^Q)  + 1 


The  flow  is  computed  at  <7.,  + 1 equally  spaced  points  on  [3  =>  1,  then 
on  p 3 2,  {3  = 3,  and  so  on. 

The  partial  derivatives  in  the  curvature  expression  eqns.  (2.7)  and 
(2  .lit-)  are  approximated  by  linear  difference  quotients.  Thus  (see  Fig.  3.U) 

(9x/0Of)p  ^ (xp , -xp) / (°rp , - 5 (du/0Q)p  ~ (up,-up)/(°<  pi-**p)  * 
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For  every  nozzle  point  and 
are  computed,  and  <T=X'-“^q  is  recorded. 


2q\p.  Xcx 


varies  mono tc-nic ally  along 


the  nozzle;  and  when  it  changes  sign  at  some  point  G (Fig.  3.1),  wc 
know  that  the  point  2,  where  the  transition  region  curvature  equals  that 
for  a simple  wave,  has  been  reached  and  passed.  Using  the  information 
at  G and  the  two  preceding  nozzle  points,  we  determine  by  three  point 
interpolation  the  variables  for<T=  0;  namely,  xE,  yE,  ru,  vE,  °<p,  . 


The  characteristic  f?=PE  (EF  in  Fig.  3.1),  on  which  x“£(c?),  (of) 

is  then  computed  in  the  same  manner  as  the  preceding  (3  = const,  character- 
istics. The  Mach  number  at  F is  the  fine]  ewit  Mach  numbers 


The  remainder  of  the  flow  is  the  cancellation  region  and  is  computed 
from  eqns . (2.9)  and  (2.10). 

In  evaluating  the  integrand  in  eqn.  (2.10) 


1 i 1 — *-c< 

L PR 

Fig.  3.5 

we  approximate  the  derivatives  . ^cx (f*)  hy  central  differences  (i.e., 
4W(«)^  r£,(cx+l)  - £ (ex  _l)J  /[2°j)  aH  points  but  the  first 

two  and  the  last  two.  For  the  first  two  points  p=[.(  ^R-$p)/(o^-oS- 


(see  Fig.  3.5),  and  for  the  last  two  points  [(^p-Sp/A^p  ”°V)  * 


r,. 


T («) 


We  approximate  the  integral  in  eqn.  (2.10),  which  we  write  as 

.r 


J (o»)  dcx  , by  Simpson's  formula  for  all  points  but- the  first- 


two  and  the  last  one:  namely, 

Cr  -1  "t 

T (c<)  ^ T (of  -2)  + 1 [_J  (of  -2)  + UJ  {&  -1)  + J (Cr)J  /3 j-  . 

T ( of  => orE ) = 0 for  the  first  point.  Tp  ^ Tp  + [(Jp  + Jp)  (Op  ~°fp)/2j 
(see  Fig.  3.5)  for  the  second  and  last  points. 

Finally,  for  the  wind  tunnel  designer’s  convenience,  the  coordinates 
of  the  nozzle  points  were  transformed  so  that  the  exit  point  F has  the 
coordinates  Xp  = 0,  Yp  = E,  an  assigned  constant.  This  calculation  places 

all  the  nozzles  on  a scale  based  on  a fixed  height  of  the  wind  tunnel  test 
section.  The  transformation  was  accomplished  by 

X - (E/V  ) (xp  - x)| 

Y - (E/V-r.)  y J 


(3.7) 
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SECTION  IV 


SUMMARY  CF  COITUTATIOTE 

Fig . U .1  indicates  the  flow  regions  and  Mach  lines  of  a typical 
nozzle  calculation. 


The  nozzles  computed  for  the  Wind  Tunnels  Branch  are  summarized 
graphically  in  Fig,  h .2 . Contours  were  computed  at  intervals  of  ,2b 
for  exit  Mach  numbers  ranging  from  1 .25  to  5 «5  • The  actual  numerical 
computations  arc  available  in  the  files  of  the  Data  Reduction  Unit  of 
the  Supersonic  Wind  Tunnels  Branch,  Exterior  Ballistics  Laboratory. 


The  designs  for  the  remaining  sections  of  the  wind  tunnel  walls 
("subsonic"  arid  "expansion"  sections  - £T  and  TB  in  Fig.  1.1)  were 
obtained  empirically  after  the  ENIAC  data  became  available . Boundary 
layer  corrections  were  included  in  all  the  jack  settings . 
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